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H. Fukaya 



1. Introduction 

JLQCD collaboration has been performing lattice QCD simulations with dynamical overlap 
fermions p]. The exact chiral symmetry realized on the lattice by the overlap fermions, not only 
enables us to analyze QCD in a theoretically clean way, but also makes the hybrid Monte Carlo 
updates fairly stable even in the vicinity of the chiral limit. In fact, we have succeeded in simulating 
QCD with the pion mass below its physical point [Q, ||, [|, ^j. 

A drawback of the use of the overlap fermion action is its high numerical cost. We therefore 
had to choose a rather small lattice size L ~ 1.8 fm, except for two main simulations with the pion 
masses around 290 MeV and 380 MeV, for which we set L ~ 2.7 fm. In particular, our simulation 
with the lightest pion mass ~ 100 MeV is inside the so-called e regime [g], where the correlation 
length of the pion exceeds the lattice size L. 

In the £ regime, the finite volume effects become large, so that the physics is largely dis- 
torted from the infinite volume limit. Such finite volume effects are mainly due to the pion's 
zero-momentum mode, and can be described well within chiral perturbation theory (ChPT) with 
the zero-mode treated non-perturbatively. We have found that the chiral perturbation formulas for 
the Dirac eigenmode density and that for the pseudoscalar two point functions [^j], nicely de- 
scribe our lattice data. We have extracted the chiral condensate [B], as well as the pion mass and 
decay constant 

A similar analysis using the chiral expansion is possible also for the three(or more)-point 
functions It is, however, expected that the contribution from the pion zero-mode become more 
complicated containing non-trivial combinations of (modified) Bessel functions, as is already seen 
in the one- and two-point functions. 

In this work on the pion form factors, we propose a different direction, or a greatly simplified 
way of analysis in the e regime of QCD. The key ideas are 

1. To insert non-zero momenta to relevant operators. 

2. To take appropriate ratios of them. 

Then, we can automatically eliminate the leading GiV) finite volume effects. Namely, we can 
extract the pion form factors without relying on non-trivial Bessel functions, as in the perturbative 
regime (p regime) Of course, this cancellation happens only at the leading order, and the next- 
to-leading order terms may contain non-trivial zero-mode contributions. But we will show here that 
the tree-level calculation already describe our lattice data for the vector form factor reasonably well, 
yielding a consistent value of the pion charge radius with the experiment. 



2. Two-point functions in the £ regime 

Before considering the three-point functions, let us begin with a simpler case with the two- 
point conelators to illustrate our new idea. For simplicity, we consider a two-flavor theory in a 
finite volume V = L 3 T, with a degenerate quark mass m. The boundary condition is set periodic in 
every direction. We denote the chiral condensate by £ and the pion decay constant by F. Note that 
including the (sea) strange quark is not difficult [11] and does not change the following results at 
the leading order of ChPT. 
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The £ expansion of ChPT is expressed by an exact group integration of the pion's zero- 
momentum mode, as well as perturbative non-zero mode fields, where their mass is treated as a 
perturbation. Namely, the theory is a hybrid system of an SU (2) matrix model (or a U(2) matrix 
model when the global topological charge of the gauge fields is fixed.), and massless fields. 



The pseudoscalar correlator separated by a four vector x = \t,x\ ,X2,x^) is, thus, given by [11] 




(P(x)P(0)) = X+Y -> +Z -> +..-, (2.1) 



where X,Y,Z- ■ • denote (dimensionful) constants, including non-trivial (Bessel) functions of rriLV 
arising from the pion zero-mode. The propagator-like forms are contributions from the non-zero 
modes. Note that the only non-zero momenta are summed. 

From the massless boson-like structure, it is not surprising to see that the zero-momentum 
correlator is a polynomial function of t. This is quite different from a conventional hyperbolic 
cosine function in the p regime. In fact, this special feature of the £ expansion, the existence of 
the constant term in particular, has been used for extracting the low-energy constants from finite 
volume lattice QCD [0, |L 0, BJ. However, in this work, we consider a different direction: how to 
reduce such peculiarity in the e regime. 

The answer is given by two steps. The first step is to insert a non-zero spatial momentum p : 

<*(.,p) S /^--<P W P(0)) = r 2£(p)sinll ' (£(p)7 . /2) c°sh(£(p)(> " T/2)) + , (2.2) 

where E(p) = |p|. Note that the constant X automatically vanishes, and the leading term has 
the same hyperbolic cosine form as in the p expansion, except for an over-all coefficient Y. The 
second step is to take the ratio of them with different momenta: Cpp(t,p) /Cpp(f,p'). Then we can 
eliminate Y which contains Bessel functions. 

Note that the sub-leading terms, expressed by ellipses, still contain the contribution from the 
zero-momentum mode. But it is suppressed typically by ~ 1 /4nF 2 \/V. 

In this way, contribution from the zero-momentum pion mode can be eliminated. It is reason- 
able that having a non-zero momentum reduces the information of the zero-mode. 

3. Three-point functions in the £ regime 

Next, we consider three point operators all put at different points : (x ^ y ^ z). In the e 
expansion, the pseudoscalar-vector(in 0-direction)-pseudoscalar function is given by 

(P(x)V (y)P(z)) =a\Y,J[ (M*-y)+jri>-*> 

1 1 (j„ n _i_ i„i \ p ip{x-y) pip 1 (y-z) 
+B V L^L { — P ° ] p2p a M(P-P') 2 ) + -, (3-D 

where Fy denotes the vector form factor, and A,B, ■ ■ ■ denote (dimensionful) constants, including 
the non-trivial contributions from the pion zero-mode. Note in this case, that the constant term does 
not exist because a "vector" cannot be made from the zero-momentum mode alone. 

In the following, we consider inserting the initial (spatial) momentum p,- to P(x), and final 
momentum p/ to P(z), splitting the discussion into three different cases. 
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Case 1 : p,-^0,p/^0 

In this case, it is straightforward to obtain 

C 3 ^ P (t,t';Pi,Pf) = [ d'xe-™* [ d 3 ze^r z (P(x)V (y)P(z)) =BF v (q 2 



x[£( Pi )+£(p/)] 



cosh(£(p,)(f-r/2)) \ / cosh(E(p f )(t' -T/2)) \ 



+ 



(3.2) 



2E(pi)smh(E(pi)T/2) J \2E(p f )sir±(E(p f )T /2) J 

where q 2 = (p; — P/) 2 — (E(pi) — E(pf)) 2 , and we have assumed t =xo — yo < T/2, and t' = 
yo-zo< T/2, and used sinh(£(p,-)(f - T/2)) ~ -cosh(£(p,-)(f - T/2)) + &(e E ^'- T / 2) ), 
and a similar relation for sinh(£"(pj)(?' — T/2)). 

Therefore, the ratio of the three-point and two two-point correlators, 

V- Cpyp{t,t'\Pi,Pf) 



R v (t,t';\pi\,\p f \,q 2 ) 



N, 



1 

3pT 



p/Up/I fixed IP'ilP/l.r 



E(pi)+E(p f ) 



N, 



2pt 

p ( .| fixed|pi 




p/ | fixedlp/ 



(3.3) 



is a good quantity to extract the form factor. Note that the rotationally symmetric correlators 
are averaged (A^", 3p ! , , and A^", 2p ! denote the numbers of correlators in the summations). 

° v lPililP/1 iPil 

Case 2 : p, ^ 0, P/ = (or Pi = 0, P/ + 0) 

Let us next consider the case with p^ = 0, but the initial momentum is kept non-zero, (the 
opposite case with p, = 0,p/ ^ 0, is obtained by simply replacing Pf = — Pi). 



In this case, the contribution from the first term of Eq. ( p.l| ) remains. It is, however, easy to 
cancel it by taking a subtraction: 

A t ,C$ P (t,t';pi,0) = C^p(t/-p h O)-CpXp(t,t ie r,Pi,0) 

cosh(J5(pj)(f-r/2)) 



BF v (q 2 



2E(p/)sinh(£(pj)r/2) 

+E{pi)T[h 



where 



V T ) 



1 



- -h 



\ T 



+ ■ 



(3.4) 



(3.5) 



and h\ (t) is its T-derivative. The reference time ? re f can be taken arbitrarily. 

For having a long-range conelation due to the zero-mode, and a periodic temporal extent, the 
correlator may have contamination from the (vac| Vo|7T7r) matrix element wrapping around 
the lattice. In fact, its contribution is not negligible at the leading order of ChPT But includ- 
ing the next-to-leading order corrections, it is suppressed by a factor ~ e - 2m *k - T / 2 \ t an d, 
therefore, neglected here. 
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Using C$(t,0) oc Th x (t/T) and its ^-derivative defined by dC PP \t,0) « h\(t/T) in the e 
regime, let us define 

£ (C$ p (t,t';pi,0) - C% P (t,t Kf ;pi,0) 



T 3pt 



i , / , i, ^Id I fixed|p,|,a 2 

4(M';|P;|,0,<7 2 ) = — — ; : : 7, (3-6) 



N, 



£ C PP \t, Pi ) -A t ,dcj$(t , ,0)+E(p i )A t ,cj$(t',0) 



2pt 

p ; | fixed|p; 



Ry(t,t ,0,0,0) = ; ^ - , (3.7) 



where A t ,dC 2 pp \t> ,0) = dC PP \t\0)-dC 2 pp \t Kh 0), and \,C PP \f,0) = C P p P (t>,0) -Cj$(fa,0), 
which is again expected to be BFy(q 2 )/Y 2 at large separations of t, and t'. Here, q 2 = p 2 . 

Case 3 : p,- = pj = 
Now it is not difficult to see 

A,A,/Cffi P (f/;0,0) 
-A t C%(t,0)A t ,dcg(t',0) - A t dC PP \t,0)A tl C Pp l (t',0) 

where we have defined A t A t <Cp§ P (t,t';0,0) = C 3 $ P (t,t' ;0,0) -C P ^ P (t,t iei ;0,0) 

—C 3 p y P (t rei -,t';0,0) +C P Y P (t Te i,t Te f,0,0), is a good quantity to extract BF v (0)/Y 2 . Fortu- 
nately, there is no contamination from the (vac|Vb|7T7r) matrix element in this correlator. 

Since Ry, R v , and R v share the same overall factor, B/Y 2 , one can eliminate it by taking their 
ratios. Noting Fy (0) = 1 , we find the following two independent extractions 

( , 2) s R v {tA\vi\y f W) R v (tA\ Pi \Aq 2 ) 
v H ' R 2 (t,t' ';0,0,0) VV ; /?|(?,f';0,0,0) ' 

with ? re f = T/4 = 12 (on our lattice), are numerically good choices, in that the signals look clean, 
while the contamination from the excited states is expected to be small. 



4. Preliminary lattice results 

Our numerical simulations are performed with the Iwasaki gauge action at j3 =2.3 and 2+1 dy- 
namical flavors of overlap quark action on a 16 3 x 48 lattice. The lattice cutoff \ja = 1.759(8)(5) GeV 
(a ~ 0.1 12(1) fm) is determined from the £2-baryon mass. The physical lattice size is L ~ 1.8 fm. 

In this work, we focus on an ensemble with the smallest up-down quark mass, ma = 0.002. 
This value roughly corresponds to 3 MeV and the pion mass is m K ~ 99 MeV ^ which is below 
the physical point. For the strange quark, we choose its mass almost at the physical value, m s a = 
0.080. Note in this set up, the pions are in the e regime (m K L ~ 0.90). 

In the Hybrid Monte Carlo (HMC) updates, the global topological charge of the gauge field 
is fixed to Q = 0. Since its effect is encoded in the pion zero-mode, the Q dependence should not 
appear in the ratios of our correlators. 

For the computation of the correlation functions, we use the smeared sources with the form 
of a single exponential function. To improve the statistical signal, the so-called all-to-all propa- 
gator technique is used: the low energy part of the correlator is separately calculated by the 160 
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Figure 1: F v (t,t',q 2 = -0.11[GeV 2 ]) (left panel) and F^tj' ,q 2 = -0.40[GeV 2 ]) (right). 
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eigenmodes of the Dirac operator and averaged over different source points, and the higher-mode 
contribution is estimated stochastically by the noise method with the dilution technique. 
Details of the numerical simulation will be reported elsewhere. 

Figure[T] shows our lattice ^ ata °^ Ey(t,t' ,q 2 ) at (p?,p 2 ,q 2 ) = (2,1,1) (left panel) and Fy (t,t',q 2 ) 
at (p?,p 2 ,q 2 ) = (1,0, 1) (right panel). Here, the momenta are denoted in units of 2k /L. To estimate 
E(p), we have used the dispersion relation E(p) = y / p 2 +~raf. Each set of data shows a plateau 
around t,t' = 5-9 from which we extract Fy(q 2 ) by a constant fit. 

The q 2 dependence of Fy{q 2 ) is presented in Fig. ^| For comparison, we have also plotted data 
in the p regime (at m = 0.015) []l2|]. Our new data in the £ regime shows a steeper slope near the 
origin, which results in a large value of the pion charge radius, 



(r 2 )v = 



dF v (q 2 
' dq 2 



0.53(4) fm 2 [preliminary]. 



This value is obtained by fitting the data to a function Fy(q 2 ) = 
(solid curve in Fig. ^) with an input of the rho-meson mass m p 



(4.1) 



+ ai q 2 +a 2 (q 2 ) 2 [0], 



l-«"7 m p 

= 770 MeV (a fit without rho- 
meson contribution (dashed curve) is also tried). This result is larger than the experimental value 
(0.452(1 l)fm 2 ) and confirms the existence of the strong (logarithmic) curvature of the pion charge 
radius near the chiral limit, as shown in Fig. ||. 



In this work, we have demonstrated how to cancel finite volume effects in the e regime. In- 
serting momenta to the operators, and taking appropriate ratios of them, we can eliminate the 
contribution from the pion zero-mode. This cancellation occurs only at the leading order, and there 
should be non-trivial corrections from NLO terms. But it is suppressed by a factor ~ l/4nF 2 \/V 
(~ 7% on our lattice). Our tree-level analysis of the vector pion form factor on the lattice in the 
e regime confirms the above observation and our preliminary result for the pion charge radius is 
consistent with the experiment, showing the existence of a diverging logarithmic curve. 

In this work, we have focused on the calculation of the vector (or electro-magnetic) form fac- 
tor. But in principle, this method can be applied for any other form factors. 
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